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SUMMARY 


Nonlinear elements are sometimes added to linear control systems 
in order to improve the response of the system to an arbitrary input. 

This can he done in different ways , one of them being the variation of 
the coefficients of the differential equation describing the system 
before the nonlinear elements are added. This variation of the coeffi- 
cients may be done in a continuous or in a discontinuous way. In the 
present paper a discontinuous variation of the coefficients is studied 
in detail and investigated for practical use. 

The nonlinear feedback is applied to a second— order system. From 
former analytical considerations the process of control is visualized 
as establishing an ensemble of linear second-order differential equations 
(some with stable and some with unstable homogeneous solutions) and 
switching from one equation to another so as to maintain small instanta- 
neous error for relatively arbitrary inputs. Physically, this control 
process is realized with a linear second-order control system to which 
have been added possible discrete combinations of proportional and deriv- 
ative feedback. The particular combination of feedback employed at any 
instant is determined by a feedback switching circuit which is in turn 
operated by sensed binary information obtained from the output, output 
derivative, error, and error derivative (namely, the signs of these vari- 
ables). Techniques that are common to the digital computer field are 
used to implement this switching circuit. 

Once physical realization is completed, simulation techniques are 
used to study and evaluate the performance of the nonlinear control system 
and to compare it with a linear system for a wide variety of inputs . In 
addition, the effects of physical imperfections that are likely to be 
encountered in any application of the control theory are considered (e.g. 
switching delays and acceleration limits). 

An analysis of the experimental results shows that this type of non- 
linear control system performs better than a linear control system having 
a natural frequency 15 times greater. For this comparison, performance 
is evaluated in terms of the average value of the magnitude of the instan- 
taneous error for band-limited inputs. Further, in contrast with the 
linear system, the nonlinear system performance is virtually independent 
of variation in the damping factor of the system. 
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A preliminary extension of this type of nonlinear control concept 
to higher order systems is presented. Experimental results are given 
for a third-order system. These results show that just as in the second- 
order case the nonlinear system performance is better than that of a 
comparable linear system. 


INTRODUCTION 


With the demand for more exacting performance, more emphasis has 
been placed on nonlinear aspects of control systems. The term ’’control 
systems” can be interpreted to include active networks and feedback 
amplifiers as well as servomechanisms. From the standpoint of analysis, 
unintentional nonlinearities have to be taken into account to explain 
performance. From the standpoint of synthesis, intentional nonlinearities 
have been introduced to improve performance. However, up to this date 
only in special cases have advancements been obtained in the field of 
non linear control systems . 

The designs of nonlinear control systems have inherent advantages. 

One advantage is that the response of a nonlinear system at a certain 
time can be made less dependent upon past response than can a linear 

system of comparable power-handling capability.^ This means that the 
nonlinear system can be made to follow more arbitrary classes of inputs 
with less dynamic error than the comparable linear system. Another 
advantage is that the mathematical difficulties encountered may actually 
be conducive to consideration of more realistic criteria of performance. 

In the nonlinear realm it is essentially as easy to invoke a criterion 
such as the minimization of instantaneous error for nonstationary random 
inputs as it is to invoke the largest possible flat amplitude response 
for sinusoidal inputs. 

In the present paper a control system of second order, which was 
first suggested by Flugge-Lotz and Wunch on the basis of analytical 
studies (ref. l), is investigated. The physical realization of this 
system and its performance are studied in great detail. 

This investigation was conducted at Stanford University under the 
sponsorship and with the financial assistance of the National Advisory 
Committee for Aeronautics. The authors wish to thank Dr. A. M. Peterson 
of the Electrical Engineering Department of Stanford University for his 
continued interest and his most helpful advice on the electronic problems 
which were encountered during this investigation. They also wish to 
thank Mr. G. S. Bahrs for his useful suggestions for a special transistor 
switching circuit. 

^In linear theory, the impulse response or the autocorrelation func- 
tion of the system gives an indication of how past response is weighted. 
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SYMBOIS 

A peak-to-peak amplitude of input 

a,b,c parameters defining a system 

a,b,c parameters in differential equation for control servo 

B viscous damping of motor and reflected load referred to motor 

shaft 

D linear damping factor 

e instantaneous error 

e a = de/dr 

H(p) transfer function, l/(ap + l) 

J inertia of motor rotor, gears, and reflected load 

j - 

^1*^2 g a in constants 

1% constant of proportionality between output velocity and back 

electromotive force 

M constant depending on initial conditions 

p operator, d/dt 

T repetition rate or period 

t time 

maximum allowable switching delay 
V input voltage 

x input into system 

y output from system 

approximation of output 
a symbol used to denote different constants 


Pm = sgn (y'e) - 2 p sgn (y’e') 


k 
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(Pm'^n) j smallest values of parameters giving good nonlinear system 


performance 
-jj 3 , gP , ^7 , g7 positive constants 

7 n = "l 7 sgn ( ye ^ ” 2 7 Sgn ^ ye '^ 

6 positive constant 

e small positive quantity 

p,v coordinates introduced in appendix D 

I\ = < 3 y/dT 


p 

radius of curvature 

0 

real part of complex frequency variable a + jo) 

T 

nondimens ional time variable normalized with respect 
ov M v t 

a 

nondime ns ional frequency, o)/o) v 

O) 

frequency 

03 v 

natural frequency of undamped linear system 

( ) 

time average 

( )' 

d( )/d T 

( )" 

d 2 ( )/dT 2 

0 

Subscripts: 

equality sign in equations which describe operations 
eqs. (2) and (3)) 

d 

ideal or desirable 

e 

error 

im 

image 

L 

limit 

lin 

linear 


m,n = 0 , 1 , 2,3 
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max 

min 

nonlin 

o 


maximum or upper bound 
minimum or lower bound 
nonlinear 
optimum 


REVIEW OF LINEAR CONTROL THEORY 


It is desirable to obtain from linear control theory some useful 
concepts that can be generalized to the nonlinear case. These concepts 
are: 

(1) Operational notation 

( a ) Trans f er func t i ons 

(b) Block-diagram representations 

(2) Control criteria 

(3) Control through parameters 

At the onset second-order systems are considered. However, there is no 
difficulty in extending these concepts to higher order systems. 


Operational Notation, Transfer Functions, and Block Diagrams 

Consider a physical process or situation in which the output is 
described in terms of the input as 

a + b ft + Cy = (l) 

dt 

where a, b, and c are constants, y = y(t) is the output, and 
x = x(t) is the input. 

Utilizing the operator p = d/dt, equation (l) may be written 

(ap 2 + bp + c)y = x (2) 

(Eq. (2) reads "(ap^ + bp + c) operating on y equals operationally x.") 
Formal solution of equation (2) for the ratio of output over input yields 
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by definition the operational transfer function for the system. Thus 

K 2 __ — i (3) 

x 2 

ap + bp + c 

The operational block diagram for the system is obtained by placing 
inside a box the transfer function, equation ( 3 ). Coming into the box 
is the input j going out of the box is the output (see sketch a) . 


X 

1 

y 

Input 

p 

ap + bp + c 

Output 


Sketch a. 

Here transfer functions and block diagrams are utilized merely as 
shorthand operational notations for differential equations . This is 
opposed to the Laplace transformation viewpoint where transfer functions 
(and thus block diagrams) have the properties of functions of the complex 
frequency variable, p = 0 + jco. The reason for stressing this inter- 
pretation is that shorthand (operational) notation has proven useful in 
the transition to nonlinear control whereas the Laplace transformation 
viewpoint (e.g., synthesis in the complex frequency plane in terms of 
poles and zeros) has not. 


Control Criteria 

To gage the performance of an actual system an ideal or desirable 
system is usually established as a straight through connection (i.e., 
y = x as denoted in sketch b, a block diagram of an ideal system). 


x 

Input 



y 


Output 


Sketch b. 


Comparison between the desired output y^ and actual output y is 
accomplished by utilizing the instantaneous error: 


e = (yd - y) = ( x - y) 


NACA TN 3826 


7 


A control criterion or criterion of performance is defined as the mini- 
mization of some property of the instantaneous error e for a given 
class of inputs. The minimizing process can be exact (i.e., resulting 
from a variational formulation of the problem) or approximate. 


Control Through Parameters 

In linear systems the process of control is usually physically 
obtained by applying feedback and/or compensation to the system that is 
to be controlled. A control criterion is realized (as closely as possi- 
ble) by adjustment of these applied quantities. The concept of control 
through parameters is an interpretation of this control process in terms 
of the differential equation describing the process. A simple example 
illustrates this concept. 

Consider the position control servo shown in figure 1. The uncon- 
trolled (open-loop) system consists of an amplifier, armature-controlled 
motor, gear train, and load. Closed-loop operation is obtained by uti- 
lizing proportional and derivative feedback. The gain constants K-^ 

and K 2 are adjustable. Armature inductance has been neglected. From 

the block diagram the differential equation for the open-loop system may 
be written: 




+ 


I B \dy = 

\ K l v dt 


x(t) 


Similarly, the closed-loop differential equation is 


/_J_\dfy + 

VWdt 2 


/ B + KgU fly 

V K lkm M 


+ (l)y = x(t) 


In either case the differential equation is of the form 


a + b + cy = x(t) 


( 4 ) 


( 5 ) 


( 6 ) 


where a set of three parameters a, b, and c completely characterizes 
the system. It is possible, then, to view the process of control in 
terms of these parameters. One starts with a parameter set (a,b,c) 
defining the uncontrolled system. A control criterion yields an optimum 
parameter set (a,b,c) Q . Control (feedback and/or compensation) is intro- 
duced ideally making it possible to adjust (a,b,c) to (a,b,c) D . In the 
above example the gain constants K]_ and Kg afford this adjustment. 
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This adjustment of the coefficients may he done in a continuous or 
a discontinuous way. In reference 2 Schmid and Triplett have described 
an interesting and efficient way to vary the coefficients of a basically 
linear system continuously. 


NONLINEAR CONTROL 


Transition To Nonlinear Control 

In the preceding section it has been mentioned that the process of 
linear control of second-order systems may be visualized as the adjust- 
ment of the parameter set (a,b,c) to the set (a,b,c) Q . The term optimum 
was used in the sense that some criterion of performance was approached 
as closely as possible. 

It seems logical in the transition to nonlinear control still to 
hypothesize control through parameters . The transition is obtained by 
allowing the parameters to become functions of the output y(t) and the 
input x(t); that is, 


a — *a(x,y) 
b — >b(x,y) 
c — >c(x,y ) 

The mathematical description of the system is now 

a(x,y) + h(x,y ) + c(x,y)y = x(t) ( 7 ) 

Without knowing the specific nature of the functions a(x,y), b(x,y), 

and c(x,y) it may be seen that equation ( 7 ) is a nonlinear, inhomoge- 
neous, and/or nonautonomous differential equation. Mathematically, 
little in general can be said about the solution of equation ( 7 ) given 
the function set a(x,y), b(x,y), and c(x,y). It seems, then, even 
more hopeless to attempt a synthesis problem which involves both finding 
the function set [a(x,y), b(x,y), c(x,y)] 0 for a specified control 

criterion and then physically realizing the system described 
mathematically . 
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Nonlinear Control Theory 

One analytical attack on the nonlinear control-system-synthesis 
problem has been made by Flugge-Lotz and Wunch (refs. 1 and 3 to 5 ). 

They suggested varying the coefficients a, b, and c, not continuously, 
but discontinuous ly. That means that for t^ < t < t 2 there is one set 

of coefficients, for t 2 < t < tj there is another set of coefficients, 

and so on. The different sets of coefficients are chosen in advance, but 
the times t^ for change from one set to another are determined by the 

value and the decrease or increase of the deviation (x - y). In other 
words , the system is linear in any interval t ^ < t < but is non- 
linear in the whole. The transition at any switching time occurs 

with continuous values of y(t) and dy/dt, but discontinuous values 
of (d 2 y/dt 2 ). 

Phase-plane techniques were used for studying appropriate sets of 
coefficients and the appropriate dependence of the switching times on the 
deviations. ^ The authors succeeded in finding a switching rule which 
assures good performance in a multitude of cases. Their control system 
is mathematically described in the following way: 

2D(i + (1 + rjy -x(t) (8) 


where 

x(t) input 

y ( t ) output 

D linear damping factor (when = y n = oj 

t nondimens ional time variable normalized with respect to co v ; 

that is, t = oo v t 


2 For details the reader is referred to references 1 and 3 to 5 . 
Reference 1 contains the ideas but is so condensed that the inquisitive 
reader will find it useful to read references 3 to 5 , of which refer- 
ence 5 is probably the most accessible. Figure 3, p. 12, and figure 30, 
p. 70 , of reference 5 will help in getting acquainted with the phase- 
plane trajectory of an output. Some of the original studies are described 
again later in the present paper when the performance of the system is 
discussed. 
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CD 


natural frequency of undamped linear system, D = 0 m = y n = 0 

Pm = -qP s s n (y’ e ) - 2P sgn (y' e ') 
m = 0, 1, 2, 3 


7 n = s S n (ye) - 2 7 sgn (ye') 


n = 0, 1, 2, 3 

-jj3 , 2 (B , -^7 , 2 7 positive constants 

/ \ f J+l for f > 0 
sgn (f) = _ = |_ x fQr f < 0 


e instantaneous error, e=^y^-y)=(x 

( )’ = 8( )/8t 


y) 


Actually, equation (8) is a normalized form of the control equation 
derived by Flugge-Lotz and Wunch. However, the notation has been some- 
what changed. See appendix A for a comparison of notations and the 
normalization involved. 

The subscript convention is 


qP + 2P 

73 = i7 + 2 / 

lP " 2P 

?2 = 17 ~ 2 / 

-qP t 2 P = "P2 

7i = -q7 + 2 7 = - 7 2 

“qP ~ 2P = “P3 

= “l? “ 2 7 = “ 7 3 


(9) 


p 0 

Properties of Equation (8) 

The following properties of equation (8) are noteworthy: 

(1) Equation (8) is a piecewise linear but overall nonlinear differ- 
ential equation. 

(2) The parameters (3^ and q n are stepwise switching functions 
of time (their implicit variable). This property is illustrated in fig- 
ure 2. 
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(3) The time of switching and the particular combination of the 

parameters and y n employed at any instant are explicit functions 

of the output y and the input x. Specifically, they are determined 
by quantized information derived from the output, output derivative, 
error, and error derivative, namely, the sign of the products sgn (y'e), 
sgn (y'e 1 ), sgn (ye), and sgn (ye'). 

(4) There are basically l6 m,n subscript combinations and thus 

l6 7 n parameter combinations . However, a detailed study of property 

(3) shows that only 8 are allowed. The allowed combinations may be 

m = n = 0, 1, 2, 3 


or 


m ^ n 
m + n = 3 

The reason for the "exclusion principle" on coefficient combinations 
stems from a desire to obtain mirrored-image outputs for mirrored-image 
inputs ; that is , 

y im ( T ) = -y( T ) 

is desired when 

x im (T) = “ x(t) 

( 5 ) The control criterion that was employed in obtaining the func- 
tional dependence of (3 m and y n was 

| yd - y| = | x - y| = |e| < e 

where e is a small positive quantity. This criterion of maintenance 
of small instantaneous error between desired output and actual output 
enabled reduction of equation (8) to the approximate autonomous differ- 
ential equation (see refs. 1 and 3 to 5 ) 

y" + 2D(l + 3 m )y ' + 7 n y = (x - y) 

« 0 ( 10 ) 

Thus phase-plane techniques could be employed to find the functional 
dependence of p m and y n . 
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(6) Once the linear damping factor D is fixed the process of con- 
trol is obtained by switching parameters p m and 7 . Equation (8) 

consists of an ensemble of eight (see property (A)) linear differential 
equations with constant coefficients. The process of control may be 
visualized as the switching from one member of the ensemble to another. 
This switching is determined by quantized information derived from the 
input and output (see property ( 3 ))- From another point of view (con- 
sistent with the approximation described in property (5))j the output y 
is to be forced to satisfy two conditions simultaneously, that is, both 
sides of equation (lO). This is approximately possible by switching to 
various and 7 n parameter combinations and can be visualized as 

the process of switching to various phase trajectories of equation (lO) 
in the phase plane of y* against y. 


Discussion of Nonlinear Control 


In the section entitled Transition to Nonlinear Control” a logical 
transition to nonlinear control systems utilizing the concept of control 
through parameters is suggested. However, mathematical difficulty hampers 
the development of this approach. In the next two sections a particular 
nonlinear control theory is presented. This theory constitutes the first 
step in the synthesis of a nonlinear control system which obtains control 
through parameters. Since the functional dependence of the parameters 
has been established, the problem is reduced to finding a set of five 
(constant) parameters (D, ^ should be appreciated, how- 


ever, that even the optimization of this five -parameter set cannot in 
general be accomplished analytically because of the overall nonlinear 
nature of the problem. 


Aside from questions on the analytical optimization of parameters 
in equation (8), there are equally important practical questions such as: 

(1) Can a useful control system that is described by the nonlinear 

differential equation (eq. (8)) be realized? 

(2) If the system is realizable, what is its physical nature? 

( 3 ) If the system is realizable, how does it compare in performance 

and complexity with a ’’good” second-order linear control system 

There are then mathematical difficulties on the one hand and physical 
difficulties on the other. The mathematical difficulties could be handled 
by numerical methods of integration of the differential equation (e.g., 
utilizing a digital computer). However, this would give little insight 
into the nature of a system that is controlled through discontinuous 
variation of the parameters and 7 n . It has been found advantageous 
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to investigate the physical questions first and then to utilize simulation 
techniques (analog computer) to investigate the analytical properties of 
this type of control. 


PHYSICAL REALIZATION 


Physical Model 

It is desirable to study the nature of a physical control system 
that is described by an ensemble of eight linear differential equations 
with control being accomplished by switching from one member of the 
ensemble to another. To do this equation (8) is rearranged as shown 
below: 

^ + 2D | + y = x < T ) - ( 2D Pm ^ + 7 n y) (id 

or in operational notation 

(p 2 + 2Dp + l)y = x - (Slipup + 7 n )y (12) 

Forgetting for the time being that the parameters and y n are 

actually functions, one interpretation of this operational equation and 
thus of equation (8) is shown in figure 3* 

Figure 3 can be modified to take into account the fact that 
and 7 n are stepwise switching functions of time, their implicit variable 
(i.e., 0 m and 7 n can each take on four discrete values). This is 

shown symbolically in figure 4. The explicit functional dependence of 
the parameters and 7 n has not yet been given and is thus indicated 

as a switching logic of undefined character. 

Utilizing the block diagram of figure 4, the physical interpretation 
of the nonlinear control system described by equation (8) is quite 
straightforward. This system consists of: 

(1) A linear feedforward portion. This portion could be a linear 
control system in itself (e.g., the simple position control servo of 
figure 1 and equation (5)). 

(2) A feedback switching circuit comprised of: 

(a) Four discrete values of proportional feedback y n 
(two positive and two negative as shown in fig. 2) 
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(b) Four discrete values of derivative feedback 2D3 m 
(two positive and two negative as shown in fig. 2) 

(c) A switching logic which at any instant determines the 
particular combination of derivative and proportional feedback 

2DP m ,7 n em P lo y ed 


Switching Logic 


In this section it is shown that digital-computer techniques can 
be utilized to establish the switching logic for the feedback switching 
circuit mentioned in the previous section. 


Recall that the parameters and y n 


tions; that is, 


have been defined as func- 


Pm = -J3 sgn (y'e) - 2 P sgn (y'e') 
m = 0, 1, 2, 3 

7 n = -p7 sgn (ye) - 2 7 sgn (ye') 
n = 0, 1, 2, 3 


(13) 


where the subscript convention has been given by equations (9). Equa- 
tions (13) determine the switching logic. Thus, for example, P^ is 

chosen when (y'e) < 0 and (y'e') < 0 and y Q is chosen when 

(ye) > 0 and (ye 1 ) > 0, so that the combination Pj7o is chosen when 

(y'e) < 0, (y'e') < 0 , (ye) > 0, and (ye') > 0. At this point it 

appears necessary to form the products ye, ye', y'e, and y'e' and 

then to find the sign of these products in order to establish the 
switching logic. Physically, however, the process of multiplication is 
to be avoided if possible. That there is a possibility of avoiding 
multiplication may be gleaned by realizing that 


gn (ab) = sgn (a) sgn (b) 


since 
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Thus equations (13) may be rewritten as 

Pm = -sgn (y')j^p sgn (e) + 2 p sgn (e'fj 
m = 0, 1, 2, 3 

7 n = -sgn (y ) |^7 sgn (e) + 2 ? sgn (e'fj 
n = 0, 1, 2, 3 

Again the subscript convention is defined by equations (9). 

Thus, utilizing equations (13), for example, p^ is chosen when 

y' > 0, e < 0, e’ < 0 

or -when 

y' < 0, e > 0, e' > 0 

7 0 is .chosen when 

y>0, e > 0, e 1 > 0 

or when 

y<0, e<0, e'<0 

so that the combination (P^^o) is chosen when 

y > 0, y' < 0, e > 0, e' > 0 

or when 

y < 0, y' > 0, e < 0, e' < 0 

From this example it can be seen that it is not necessary to find the 

signs of products but rather that it is sufficient to find separately 
the signs of y, y', e, and e'. 

Since the sign of a variable is quantized binary information of the 
variable, it is convenient to utilize digital-computer techniques to 
further the switching logic. This may be done as follows: 

Let the convention be adopted that y > 0 be represented by 0 
(binary zero), y < 0 be represented by 1 (binary one), and similarly 
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for y', e, and e'. If the ordered sequence is now established as 

(y, y' , e, e') 

four-digit binary logic may be employed to encode equation (l4). In 
particular, a binary coded decimal may be used (see table l). 

It was mentioned in property (4) of the section Properties of 
Equation (8)" that not all of the l6 possible p^,7 n parameter combina- 
tions were allowed under their definition. This was termed an exclusion 
principle on the allowed coefficient combinations. One of the advantages 
of the suggested binary coding scheme of table 1 is that this exclusion 
principle is built into the code. To understand this, consider the 
example of the combination 0 5 7 O given previously. In the code language 

P^ 7 q is chosen when the binary number 0100 or 1011 occurs. What is 

implied by this example is that a binary number and its complement must 
be identical (i.e., 0000 = 1111 , 0111 = 1000 ) as far as the switching 
logic is concerned. Thus out of the 16 possible four -digit binary numbers 
only the first 8 are unique. That is, in counting from 0 to 7 in a binary 
coded decimal, if complements are included then so are the other 8 possi- 
bilities, 8 to 15 (see table l). 

The allowed 1^,7 parameter combinations along with the encoded 

logic of table 1 are summarized in matrix form in table 2. Examples are 
given to illustrate the meaning of the table. In general the allowed 
subscript combinations are 

m = n = 0, 1, 2, 3 
or 

m ^ n 

m + n = 3 


Realization Completed 

Now that equation (l4) has been successfully interpreted (encoded) 
in binary-logic form (table 2), the realization of a feedback switching 
circuit utilizing this encoded logic is a typical digital-computer 
switching-circuit problem. As is characteristic of any synthesis process 
there will, in general, be many ways to design this feedback switching 
circuit. The block diagram of figure 5 shows one design that completes 
the physical interpretation of equation (8) along the lines started in 
figures 3 and 4 . 


5=0 

T 
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In figure 5 the feedback switching circuit consists of: 

(1) The four discrete values of both derivative and proportional 
feedback SD^, 7 n 

(2) A relay switching circuit that connects the proper feedback 
combination 

(3) Zero-coincidence detectors CD that drive the banks of relays 
to one position or the other depending upon the signs of the sensed 
variables 

It should be noted that, depending upon the application, other forms 
of sign-sensing devices and other switching devices such as diodes, tran- 
sistors, electronic switches, and/or magnetic amplifiers could be employed 
to obtain other realizations of equation (8). In any case the following 
properties are basic to any realization: 

(1) The signs of the four variables y, y', e, and e' are 
sensed. This may be thought of as the process of "reading in" the four- 
digit binary logic of table 1. 

(2) On the basis of the 2^ possible binary decisions the required 
feedback combinations 20^, 7 n as defined in table 2 are connected around 
a linear second-order member. 

It is important to stress that the only type of nonlinear operations 
required in the realization of this nonlinear control system are switching- 
type operations. In addition, all the switching is to be performed in 
feedback paths, which means that the switching can be done at low elec- 
tronic power levels. These practical features are definite design advan- 
tages. Thus, in summary, it can be said that this type of nonlinear 
control system is not only physically realizable but also practical from 
an instrumentation standpoint. 


EXPERIMENTAL VERIFICATION OF NONLINEAR CONTROL THEORY 


Discussion of Simulation Techniques 

Simulation techniques were chosen as an experimental mode of inves- 
tigation of performance of the nonlinear control system realized from 
equation (8). These techniques offer the following advantages: 
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(1) Proximity to the actual control system. This means that the 
same practical features with regard to instrumentation (see the section 
entitled "Realization Completed") are exploited to fullest advantage. 
Thus . just as in the actual system, the only nonlinear device required 
for the simulated model is a binary logic feedback switching circuit 
(see appendix B for details). The linear portion of the system is simu- 
lated on an analog computer. Here the only operations required are sum- 
mations, two integrations, and one differentiation. These are all 
operations which an analog computer does well. It can be said then that 
the accuracy to which the simulated model simulates equation (8) depends 
primarily upon the realized feedback switching circuit. The most essen- 
tial type of imperfection to be expected in this switching circuit is 
time delay in switching. Exactly the same type of imperfection will be 
met in the physical control system. Thus there will be more nearly a 
one-to-one correspondence between the simulated model and the actual 
system than between either and equation (8). 

(2) Convenience in experimental investigation. In order to charac- 
terize the output y of the nonlinear system completely, a set of five 
parameters 


(3, pPj 2 ^’ ~\ 7 ’ 2^) 

and the input x must be specified. In the performance evaluation of 
the system it is necessary to be able to vary these characterizing 
quantities conveniently. Simulation techniques allow this. 


Presentation of Experimental Results 

Figures 6 to 17 present experimental results obtained from the simu- 
lation studies of equation (8). Briefly, the results are presented as 
follows: 

Figures 6 to 12 compare the responses (output y and error e) of 
the nonlinear system with that of a linear system for various classes of 
inputs x. (in comparing the linear and nonlinear responses it will be 
noted that there is not exact synchronism of events because, with the 
available experimental facilities, it was necessary to obtain the two 
responses separately.) The linear system utilized is that which consti- 
tutes the feedforward member of the nonlinear system corresponding to 
the case where 3m = 7n = °* 11516 nonlinear system for figures 6 to 15 is 

3 3 = - 3 0 = 2 

32 = “3i = 3-5 
73 = -7 0 = 2 
7 2 = “7]_ = °*5 


NACA TN 3826 


19 


Figure 6 compares the system responses to sinusoidal inputs and 
f igure 7 shows the responses to triangular-wave inputs . A partially 
integrated square wave x(t) is defined as the output of a first-order 
linear system characterized by the transfer function H(p) = l/(ap + l), 
when the input (t ) is a square wave. The responses of linear and 

nonlinear systems to this type of input are given in figure 8. The 
responses of the systems to small sinusoidal inputs are given in fig- 
ure 9* Figure 10 gives the responses to clipped sinusoidal inputs. Fig- 
ure 11 shows the response of the nonlinear system to sinusoidal inputs 
that have been displaced with a direct-current component. Figure 12 
shows the responses of linear and nonlinear systems to a triangular -wave 
input whose periods and amplitudes are randomly modulated. 

Figures 13 to 15 deal with the effects of imperfections that are 
likely to be encountered in the actual control system. Figure 13(a) gives 
the results and data of an experimental investigation on the effects of 
switching delays due to threshold in sensing the sign of the error, 
sgn (e) foi a triangular -wave input. The experimental results for a 
constant 9-volt input (see fig. 13(b)) are given below: 

Threshold, mv 14 26 36 44 

Peak-to-peak error, mv 44 124 220 290 

Figure 14 shows the effects of placing progressively smaller limits on 
the acceleration of the nonlinear control system. For each value of y" 
limit considered, the output, output derivative, output acceleration, 
and instantaneous error are shown. The effects of a velocity limit on 
performance of the nonlinear and linear systems are compared in figure 15- 

Figure l6 presents the responses of some special cases of the non- 
linear system (special with respect to the choice of the 3 m and 7 n 

parameter values ) . In this figure the system responses to a triangular- 
wave input are given for four different parameter sets . The parameters 
pertinent to these results are listed in table 3 for easy reference. 

Figure 17 gives the response of the nonlinear system possessing a 
low linear damping factor D = 0.1. The responses are for a triangular- 
wave input whose frequency was varied in the same manner as that of fig- 
ures 7(a) and 7 (c). 

Detailed discussions of these results are given in the section 
entitled "Discussion of Results." 
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DISCUSSION OF RESULTS 


Performance Evaluation of Sinusoidal Inputs 

To complete the synthesis of the nonlinear control system that has 
heen derived from equation (8), it is necessary to choose the magnitudes 
of the parameters that characterize the system, that is, 

or (D,Pj,P2>73j 72) ' ^ is ex P ec ' te< ^ that the performance of the system 

depends on the choice of these parameters. For studying their influence, 
simulation techniques proved to be very convenient. Experimental results 
were given in figures 6 to 12 where the response (output y and error e) 
of the nonlinear system was compared with the response of a linear system 
for a variety of inputs. This gives the possibility of establishing the 
properties of the nonlinear system not only by itself but also with respect 
to a linear standard. The linear system employed for this purpose was 
that which constitutes the feedforward member of the nonlinear system 
(i.e., ^ = 7 n = °) • 

These experiments allow parameter values for good performance of the 
nonlinear system to be found. Analytical and practical considerations 
that aid in the optimization are treated later in the section entitled 
"Choice of Parameter Values." 

The sinusoidal responses of the nonlinear and linear systems are 
compared in figure 6. Here, the frequency range considered was 
0.1 ^ ft = ao/aXy = 2. The peak-to-peak input amplitude was 20 volts. These 

results show that the nonlinear system reproduced the sinusoidal inputs 
up to the frequency ft = 1.4 with virtually no instantaneous error when 
compared with that of the linear system. For higher frequencies the error 
for the nonlinear system increased rapidly to the same order of magnitude 
as that of the linear system. 

Substantially the same comparative performance was displayed by the 
two systems when the input was a smaller (4-volt peak-to-peak) sinusoid 
as is shown in figure 9* 

Figure 11 gives the response of the nonlinear system to a 20-volt 
peak-to-peak osculating sinusoid over the frequency range 0.1 ^ ft S 0.8. 
This is a severe type of input for the nonlinear system since both x 
and x ! simultaneously go to zero. This implies that y and y f are 
also small so that in equation (8) the discontinuous variations of the 
parameters cannot be so effective in determining the acceleration 

y", since 


y = 


x - 


[pC 1 + fta)y' + i 1 + 7n)y] 


(15) 


NACA TN 3826 


21 


From the figures it is seen that the nonlinear system did have some dif- 
ficulties near the osculating regions; further, the peak error increased 
as the frequency was increased (see appendix C). As might be expected 
from equation (15) this error for small values of y and y' can be 
reduced by increasing the magnitudes of the parameters (this will be 
discussed in more detail in the section entitled "Choice of Parameters"). 
At any rate, by comparison, the error for the present system is always 
less than that of the linear system of figure 6. (in the linear realm, 
error is independent of a shift in the direct-current level of the input.) 

Taken collectively, the 20-volt, the 4-volt, and the osculating 
20-volt sinusoidal inputs tend to form a more realistic appraisal of the 
nonlinear system performance than a single input amplitude. There is 
still much that can be learned from a detailed study of these three 
responses, but first it is desirable to obtain some sort of a quantitative 
comparison between the nonlinear and the linear system performance. One 
such comparison can be made as follows: 

(1) Assume that the input to both systems x(r) is, and has been 
for a long time, a 4-volt peak-to-peak sinusoid of frequency 

ft = cd /< n v = 0.1 

so that as far as the linear system is concerned this is a steady-state 
alternating-current input. 

(2) Determine how much the band width or the natural frequency o> v 

of the linear system must be increased in order that the time average of 
the magnitude of the instantaneous error 



| e |dT 


for the linear system be reduced to that value given by the nonlinear 
system. 

Here it is easy to show that for the low-frequency steady-state 
alternating -current case 

(16) 


llin 




max 


For an input x = |x| max e^ T the steady-state error is given by 


e = x - y = I x | max ( 1 - 


JGt 


1 + jft2D - ft c 
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For small values of 0 . this equation yields 

e - |x,_[l - (1 - J02D + . . . )] e JIlT 
- UI„ lax (3fl2D)eJ nT 


or 


but 


I e I max ~ 2Dft|x | 


max 


Him = J l e lmax sin *** = §l e| 


max 


therefore 


Ie| 


lin 


T fi|x| ni8a 


(IT) 


For the given system and input 



0.15 


For the nonlinear system and the same input (approximate calculation 
from fig. 9(a)) 


Hnonlin 0.01 


Thus, 


^ e ^ nonlin ~ 15 ^ 6 ^ Tin 

Now from equation (17), since Q = co/cn v , it is seen that in order to 

reduce |el lln t>y this factor of 15 that a>y must be increased by the 
same factor. 

In the example of the linear control system of figure 1 and equa- 
tion (5) this increase in o) v by a factor of 15 would mean that both 
the inner and outer loop gains must be increased by the order of magni- 
tude of (15) 2 since 
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K A . 

J 


(av ) 2 


This is true assuming that the damping factor remains constant. Such 
an increase in the loop gains is frequently not at all physically possible. 

Up to this point little attention has been given to the detailed 
nature of the nonlinear system response. Closer inspection of, for 
example, figures 6(a) and 6(c) shows that the output y is a function 
that links serpentine fashion (oscillates) at a very high rate about the 
input x, but still the magnitude of the error is small. In fact, it 
is necessary to inspect the error at a scale 20 times larger than that 
of the output even to notice this phenomenon. Mathematically this means 
that the functions x and y approach one another closely but that 
their derivatives differ appreciably. Physically, however, this is not 
at all undesirable as long as the magnitude of the error is small. 

(Actually for mechanical systems this property would be useful in pre- 
venting static friction. ) This fine-grained oscillating character of y 
is the very essence of the nonlinear control theory. Every time the 
error or error derivative goes through zero the parameter set 7 n of 

equation (8) changes discretely as defined in equations (13) or by the 
binary logic of table 2. The discrete changes in the parameters cause 
discontinuities in the second derivative y", which when integrated twice 
give y its serpentine character. To illustrate this point, the sketch 
of figure l8 shows samples of the superposed input and output of the 
nonlinear system. The input in this case could be that of figure 6(a) 
or 9(a). 

Circle (a) in figure 18 is a typical cycle of the oscillating char- 
acter of the output y. Commencing at Tq, the error changes sign at 

Tq, the binary number 0011 is "read into" the feedback switching circuit 

of figure and using the notation of table 2 the parameter combination 
|3 j 7 j is switched into the circuit. This causes an immediate reversal 

in the sign of the output acceleration y", so that at time the 

error derivative changes sign, the binary number 0010 is read in, and 
the parameter combination (3 2 ?2 is switched into the circuit; acceler- 
ation is still in the same direction but weaker. At time the error 

again changes sign, the binary number 0000 occurs, and the combination 
Pq7q is switched into the circuit; acceleration is in the opposite 

direction. At time the error derivative changes sign, the binary 

number 0001 occurs, and 3^1 switched in; this reduces the acceler- 
ation until at time the error again changes sign, 0011 occurs, 3^7^ 

is again switched in; the cycle is complete. Although it was not 
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mentioned at the time, figure 2 shows this sequence of combina- 
tions. The nature of the acceleration resulting from the switching can 
be seen in figure 14(a). 

The comparable switching cycle for y > 0 and y' < 0 as shown 
in circle (b) of figure l8 is: 


Logic 

0111 

0110 

0100 

0101 

0111 

Similarly, for y < 0 and y' < 0 


Parameters 

P 0 7 3 

P l 7 2 

p 3 7 0 

P 2 7l 

Po^3 

as shown in circle (c) of figure l8: 


Logic 

1100 

1101 

1111 

1110 

1100 


Parameters 

P 5 7 5 

P 2 7 2 

P 0 7 0 

Pl 7 l 

p 5 7 3 


This is the mirror-image switching cycle for circle (a) (see property (4) 
in the section "Properties of Equation (8)" and also the discussion in 
the section "Switching Logic"). 


Finally, for y < 0 


and y' 

>0 as in c 

Logic 

Parameters 

1000 

p 0 7 3 

1001 

Pi 7 2 

1011 

p 3 7 o 

1010 

P 2 7 i 

1000 

p 0 7 3 


as in circle (d) of figure 18 : 


This is the mirror-image switching cycle for circle (b). 

With this insight into the detailed behavior of the nonlinear system 
more information can be obtained from the experimental sinusoidal responses 
of figures 6, 9> and 11 that have up until now been treated from a macro- 
scopic rather than microscopic viewpoint. Along these lines, the following 
experimentally observed facts are noteworthy: 
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(1) In comparing the errors for the 4-volt and the 20-volt peak- 
to-peak sinusoidal inputs in the frequency range of good reproduction 
0.1 ^ ft ^ 0.8, it is seen that: 

(a) The magnitudes of the errors are nearly the same (see 
figs. 6(a) and 6(e) and 9(a) and 9(c). 

(b) The period of the error is generally smaller in the larger 
amplitude case. This is even more pronounced if the 4-volt peak-to-peak 
case is compared with figure 11 for y in the region of -20 volts. 

(2) Good reproduction is characterized by many e and e' switchings 
(see fig. 18) per cycle of the input. As input frequency is increased 

and the upper limit of small error is reached (figs. 6(e) and 9(e)) the 
e and e' switchings become more infrequent until there are finally 
only two of each per cycle of the input. 

The fact that the period of the error is smaller for larger inputs 
can be gleaned from equation (8) when it is rearranged as 

y" = -[2D(1 + (3 m )y' + 7 n y] + (x - y) 

or roughly 

y" « -[2D(l + Pnjx' + 7 n xJ (18) 

For a system with specified parameter values, the larger input 

and input derivatives will give stronger discontinuities in y" as the 
&m> 7 n combinations change. This implies that 7 n will change more 

often making the period of the error smaller. See figure 14(a) for the 
justification of the approximations in equation (l8) since y, y', y", 

and x are shown in this figure. (The input x was sketched in by hand.) 

In the frequency range of good reproduction the reason that the 
magnitude of the error is virtually independent of the input amplitude 
cannot be explained from equation (8) since this fact is intimately tied 
in with the imperfections in the feedback switching circuit (see section 
entitled "Effects of Switching Imperfections" for details). As measured 
from the experimental sinusoidal responses, |e| max is of the order of 

15 to 20 millivolts. It might be noted that earlier in this section the 
smaller 4-volt peak-to-peak input sinusoid was employed in the compari- 
son when it was determined that o) v should be increased by a factor 

of 15 in order to obtain the same This choice of the smaller 

input was decidedly in favor of the linear system since the linear error 
increases linearly with input amplitude. Thus, if the 20-volt peak-to- 
peak input were used, an increase in ou v by a factor of 5 X 15 = 75 would 
be required to obtain the same l e l]_i n - 
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The fact that there is a correlation between good reproduction 
(small error) and the existence of many error and error-derivative 
switchings is very important since it is the key to understanding the 
upper limits of good performance of the nonlinear system. One approxi- 
mate way to investigate this matter is to establish a deficiency between 
the output acceleration required for good reproduction (many switchings ) 
and the available output acceleration. Equation (19) gives an approxi- 
mation to the output acceleration when the system is functioning well: 

y" « -[2D(l + Pm) x ' + 7 n x ] U9) 


As an example of this approach consider the input x to be a sinusoid 
of frequency ft (fig. 6 or 9)> Since 


x 


1 max 


= ft 2 |x| 


max 


it is to be expected that the nonlinear system will have the greatest 
difficulty in the vicinity of | x | max where, from equation ( 19 ), 

I y I max ~ I ^n x I max 


Thus, if in this region the output y is going to interweave the 
input x (as is characteristic of good reproduction), then 

I v" I > lx" I 
I ■ y I max 1 A 1 max 


This inequality then places an upper limit on good performance of the 
system in response to sinusoids 

^ < ^ l^nlmax 

For the system used in obtaining the experimental results of figures 6 
3X111 |^n|max = 2 » therefore, 


a < {2 

Inspection of both figures 6(e) and 9(e) tends to substantiate the above 
result. For example, in figure 9(e) it is seen for ft = 1.2 (after the 
transient caused by turning on the input has been absorbed) that as the 
input goes through its maximum, frequent e and e' switchings stop 
and do not occur again as the frequency is increased. 


NACA TN 3826 


27 


Performance Evaluation of Inputs Other than Sinusoidal 

In the preceding section only the sinusoidal response of the non- 
linear control system was discussed in studying performance. However, 
it is easy to see that the nonlinear system will cope with any input in 
the same manner as it does with sinusoids. Thus, as long as the magni- 
tude of input acceleration does not continuously exceed the available 
magnitude of the output acceleration the switching process will commence 
and excellent reproduction will result. This type includes inputs with 
discontinuous derivatives and discontinuous inputs. The experimental 
sinusoidal responses of figures 6, 9> and 11 themselves give some indi- 
cation of these facts. For example, in figure 6(a) there was an initial 
discontinuity in the input and there were discontinuities in the input 
derivative when the frequency was changed. The results in figures 7, 8, 
10, and 12 prove further that the nonlinear system response is not 
dependent upon any specific type of input. Given in these figures are 
triangular-wave, partially integrated square-wave (including square-wave), 

clipped sinusoidal, and random^ inputs, respectively. 

It should be noted that in the literature (refs. 1 and 3 to 5) a 
marginal- type input that would present a case of indecision to the 
switching circuit of the nonlinear system is discussed. This case has 
never been encountered experimentally even when the attempt was to pro- 
duce this case. Thus the marginal- type input is not considered practi- 
cally important. 


Use of Phase-Plane Methods to Study Performance 

The phase-plane methods that were used in the original analytical 
development of equation (8) (refs. 1 and 3 to 5) can also be gainfully 

^In the strict statistical sense the probability that the inputs 
shown in figures 12(a) and 12(b) are not samples of random stationary 
time series is admittedly high because of existing experimental facil- 
ities. These inputs were obtained by random manual modulation (both fre- 
quency and amplitude) of a triangular -wave input of peak-to-peak amplitude 
A and period T where 


0 $ A ^ 20 volts 
0.1 £ = 2n/T ^ 0.8 

Random manual modulation means that the operator varied by hand both the 
frequency and amplitude controls of the input generator as randomly as 
possible. In the present investigation the comparative results of the 
linear and nonlinear responses to what appear to be band-limited random 
inputs are felt to be more important than the exact statistical proper- 
ties of the inputs . 
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employed in studying "the performance of the nonlinear control system 
derived from equation (8) once the p m and 7 n parameter values have 

been specified, (in this section it is important to distinguish between 
Pjjj. and 7 n parameter values, i.e., (3^ = 2, f$2 = etc., and 1 Vn 

parameter combinations, i.e., p^7o> etc.) 

If in equation (8) the error e = (x - y) is assumed small, then 
the output y can be approximated by portions of curves that satisfy 
the autonomous differential equation 

2^ + 2D(1 + Pn,) f£ + V ' 0 < 20 > 

dr 

where D, p m , and 7 n are defined in the section "Nonlinear Control 
Theory." Further, if dy/dr = f-,, these approximating curves are defined 
in the f^y phase plane by 

(yx 2 " Ii)* 2 = M (^i - hf 1 ^ 

where M is a constant depending on initial conditions and 

*1,2 * -Dp + Pm) t |d 2 (i + p m ) 2 - 7 n 

Equation (2l) comes from integration of the firrt-order differential 
equation 


<Hi = -2D(1 + Pm)^ - y n y 
dy Ip 

= -2D(1 + Pm ) -7nf < 22 ) 


Since the particular P m ^7 n combination employed in each point of 
the phase plane depends upon y, Ip e, and e', four approximating 

curves go through each point (see appendix D). The tangents to these 
curves (eq. (22)) indicate four directions which lie in an angular sector 
(see fig. 19). This angular sector is defined by the two extreme direc- 
tions which apply if e and e' have the same sign. The two inner 
direction values apply if e and e' have opposite signs. By super- 
posing the input x'x phase plane on this l-^y plane it can be stated 
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that at any point the tangent to the input phase curve must he included 
within the shaded angular sector (discussed above) if small error is to 
be obtained. Thus, equation (22) can be used to study performance. The 
width of the angular sector changes with the ratio y/f^ for given P m 

and 7 n parameter values. It is largest (l80°) for |q — >0 and smallest 
for y = 0 if (l + Pm) and 7 n have the same sign. However, the width 
for £q — ^0 is not significant since all phase curves (including the 
input) have infinite slope there. Thus, along the line = 0 the 

curvature of the input phase curve must be used to determine limits on 
good performance. The radius of curvature of the input must be smaller 
than that of the flattest approximating curve. 


In order to illustrate the use of these phase-plane methods in 
predicting limits of good performance consider the example of a sinus- 
oidal input x = sin At. In this case the input is represented by an 
ellippe in the phase plane and 

?-') = _fi 2 2L 
dx x' 


In the superposed xx' and y|q phase planes of figure 20 are shown 

three sinusoidal-input phase curves (i.e., three different frequencies) 
and the families of phase trajectories of equation (22) for the |Bm, 7 n 

combinations where e and e' have the same sign, (in order to avoid 
extensive computations, the parameter values indicated in figure 20 are 
those of an earlier investigation (see ref. l).) Thus, the tangents to 
two intersecting phase curves define the angular sector at that point, 
as has been discussed. By tracing the inputs with different values 
of ft it can be understood that good performance for the presented 
system can be obtained only for ft < ft Q where ft 0 is the parameter 
belonging to that^e 1 lip s e which has the same radius of curvature as the 
curve through (y,£q) = (l.O). (Note that the radius of curvature of the 

approximating curves jumps at = 0 and is smaller for ^ < 0. ) In 
the first quadrant (i.e., approaching fq = 0 f rom Tq > 0) this radius 
of curvature is given by 


" |p| r ' v 


Since at this point the ellipse radius of curvature is ft 2 /l. 


< 7 n y 
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Therefore, for the first quadrant 

n ( 23 ) 

Here it should he noted that this is the same result as that obtained 
in the section entitled ’’Performance Evaluation of Sinusoidal Inputs” 
even though the parameter va lues are different. For the system of fig- 
ure 20 7^ = 0.5; thus ft < )/0.5. 

Another example is given by the input x = 1 - e~ aT with 
d(x ! )/dx = -a. Since the smallest angular sector is at y = 0, 
a = 2D(l + P max ) determines the limits on performance. This means that 

for the system represented in figure 20 good control can be expected for 
a value of a slightly smaller than O. 75 . 

A step input is represented by x = 1 - e~ aT with a — > 00 . The 
picture in the phase plane is x = 1 - (l/cOx 1 with a — >00. For very 
large values of a this is a straight line which forms a small angle 
with x ! =0. For a — the curve degenerates to a point (x = 1, 
x 1 =0). There is no doubt that a perfect followup of a step is not 
possible because for practical reasons the line for large values of a 
does not lie in the allowed angular sector at any point it is passing 
through. 

Related to the step input is the square-wave input. A portion of 
the output phase trajectory for a square-wave input to the system of 
figure 8(e) is shown in figure 21. Figure 21 is computed with the help 
of the differential equation 


_ _ h 

dll 2D(1 + + (1 + 7n)y - x 

This equation is obtained from equation (8) by replacing d y/d t by 
£l(d£i/dy). In this example one cannot immediately use approximating 

curves for designing the output because the error is too big at the start 
of the motion. The computed diagram is in good agreement with the test 
run shown in figure 8(e). 


Effects of Switching Imperfections 

In equation (8) it is assumed that the parameters p m and y n 

change upon exact zero coincidence of any one or more of the variables 
y, y', e, and e 1 as defined in equation ( 13 ) or by the logic of 
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table 2. Physical imperfections, however, preclude this possibility. 
Thus, in the simulated model of equation (8), exact zero coincidence 
cannot be detected because of threshold effects, and relays are subject 
to time lag (both mechanical and electrical), dead zone, and chatter. 

In the present section an attempt is made to evaluate the effects 
that these switching imperfections had on the experimental result in 
order to obtain some practical design criteria for specifying switching 
requirements for good performance in other applications of this method 
of nonlinear control. 

The experimental results of figures 6, 9, and 11 can be used to 
demonstrate that relay imperfections were not important in the simulated 
model. Comparison of the sinusoidal responses in these figures has shown 
that the period of the error became smaller as the magnitude of the input 
amplitude was increased. (This was discussed in the section "Performance 
Evaluation of Sinusoidal Inputs.") From figure 11(a) in the region of 
y » 20 volts (the largest magnitude considered in all the experimental 
studies) the period of the error T e is measured as 

T e = o) v t e » 0.2 

There are four parameter switchings per error cycle (see fig. 18) . 
Assuming these to be approximately uniformly spaced, the mi n i m um time 
between parameter switchings is approximately T e /4 - Now assuming that 

the relays must be capable of closure in at least 1/5 of this time, the 
maximum allowable (real time) switching delay t^ is 

“vW ~ T e/ 20 = 10 " 2 (23) 

In the simulated model a> v = 1 radian per second. Therefore, 

tfflfl ~ 10 millisec 

As given in table 4, the relays employed were capable of closure in 
3 milliseconds or less so that they were entirely adequate for the experi- 
mental studies. 


Bie ability to disregard the relays in the evaluation of the effects 
of switching imperfection on performance leaves only threshold effect in 
sensing the sign of the variables y, y', e, and e' to be considered. 
As has been observed and discussed in the section "Performance Evaluation 
of Sinusoidal Inputs the amplitude of the error for the nonlinear system 
was relatively independent of the magnitude of the input (i.e., when the 
system is operating in the rapid e to e' switching sequence so that 
this is the minimum-error case). This constancy of the lower limit on 
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the error magnitude is caused by switching imperfections and thus thresh- 
old in sensing the signs of variables. Since under a normal switching 
sequence there are many more e and e' switchings than y or y' 
switchings and since the error is of an order of magnitude less than the 
error derivative , the primary cause of the lower limit on error is local- 
ized as threshold in sensing the sign of the error. Figure 13 shows the 
results of an investigation of this threshold effect. In all system- 
response figures except figure 13 the peak-to-peak threshold was approxi- 
mately lA millivolts. In figure 13(a) the peak-to-peak threshold was 
hb millivolts so that figures 7 (a) and 13 give a good comparison of the 
effects of these two threshold values. 


Choice of Parameter Values 

The performance of a completely specified nonlinear system has been 
discussed. The parameter values for this system, that is, 

D = 0.6 

P3 = “Po = 2 

P 2 = -Px = 0.5 

73 = -7 0 = 2 
7 2 = -7 l = 0-5 

were initially chosen in the following manner: D was first selected 

to give good linear system performance; then the smallest values of the 
3 and 7 n parameters giving good nonlinear system performance were 

chosen experimentally from a systematic variation of parameters utilizing 
the simulated system. This particular set of 3m and 7 n parameter 

values can thus be denoted as (Pm^7n)min since they establish the lower 
bound on parameter values for good nonlinear system performance. The 
physical significance of (Pm>7n) m in is that loop gains and acceleration 

requirements of the linear member are minimized since 3 m and y n are 
feedback gain constants . 

From the discussion in the section "Performance Evaluation of 
Sinusoidal Inputs" centering about equation (l9) or from the phase-plane 
methods of the preceding section it is to be expected that a general 
increase in the parameter values over (Pm> 7 n) m i n will result in improved 

nonlinear system performance by increasing the available acceleration of 
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the system or increasing the angular sectors in the phase plane. Inspec- 
tion of figures 16 (a), system 1, and 16 (c), system shows this to he 
true. In figure 16(a) the parameter values were ( 0m>7n)m-tn 85 given 

above, and in figure l6(c), system they were 


0 3 = - 3 0 = 10 
0 2 = -01 = 1 



On comparing the performance of these two systems it is noted that the 
corners of the input triangular wave are reproduced with less error by 
the system of figure l6(c). Thus far then it would seem that there is 
no upper bound on the parameter values; that is 


(<Wn) max ->~ 

Invariably, however, there will be upper bounds on the values of the 
parameters because of acceleration limits in the physical system. Fig- 
ure 14 shows the effects of placing progressively smaller acceleration 
limits on the system. Here it is seen that for |y ,l |j j > 0-3|y"| mflv per- 
formance is not appreciably affected but for values less than this good 
performance is no longer obtained so that acceleration limits definitely 
tend to determine (0m,7 n ) max - In general, then, there will be a whole 

range of values of p m and 7 n ; that is, 

(0m^n) m ax > (Pm>7 n ) > (0m->7n) m i n 

for which good nonlinear system performance results. The final choice 
must depend upon the particular application and can easily be found 
experimentally . 

There are certain special cases of the p m and y n parameter values 

that lead to simplified feedback switching circuits and thus lead to non- 
linear systems that are simpler to realize. In table 3 three of these 
are denoted as (l) no derivative feedback, (2) no e' sensing, and 
(3) no e sensing in derivative feedback loop and no e' sensing in 
proportional feedback loop. 

By making inoperative the appropriate relays in figure 5 the sim- 
plified switching circuits for these cases are easily obtained. It is 
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desirable then to inspect the performance of these special cases (in 
comparison with that of the complete system) to ascertain whether or not 
as good performance can be obtained with less complexity. Figure 16 
shows the response of these cases in comparison with that of the complete 
system. Figures 16(a) and 16(b) differ from figures l6(c) and l6(d) in 
that the general magnitude of the parameters was increased in figures l6(c) 
and l6(d). 

These results show that case (2) is not worthy of much consideration 
since the magnitude of the oscillating error is large. Cases (l) and (3), 
however , should he considered for certain applications. For example, if 
it were known that the amplitude distribution of the input was relatively 
void near zero, then case (l) would serve as well as the case of the 
complete system. Case (3) shows nearly constant percentage error so that 
it could be useful in cases where accuracy was not so important as economy 
in components . 

To this point the parameter D has received little attention mainly 
because its value (within limits) is not particularly important. It has 
been observed experimentally that D may be anywhere in the range 
0 < D < 1 even for (3 m ,7n) m i n and. performance of the nonlinear system 

is not affected. Inspection of the block diagram of the nonlinear system 

(fig. 4 or 5) shows that the physical significance of the variation 

of D for given values of P>ni,7 n is that the damping factor of the linear 

member and the derivative feedback around the linear member change syn- 
chronously. A case of more practical importance such as might arise in 
aerodynamic applications of this type of control system is the variation 
of D for given values of Dp ni ,7 n . That is, the damping factor of the 

linear member alone varies while the feedback values around this member 
remain unaffected. Even in this case it has been found that the nonlinear 
system performs well. For example, figure 17 shows the response of the 
nonlinear system 


D = 0.1 

= -3q = 12 

p 2 = -p 1 = 3 
73 = - 7 q = 12 

7 2 = -7 0 = o-5 


to triangular-wave inputs, while figure 7 shows the response of the system 
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D = 0.6 

% = “Po = 2 

02 = “0q = 0*5 
7 5 = -7; L = 2 
7 3 = -7 0 = 0.5 

to the same inputs. Comparing these results it is seen that performance 
is nearly the same for both systems. 


PRELIMINARY EXTENSION TO HIGHER ORDER SYSTEMS 


From a practical standpoint limitations in the applicability of the 
nonlinear control system described by equation ( 8 ) do not stem from the 
inability to realize the feedback switching circuit but rather from 
assuming the linear member to be of second order. In many cases a more 
realistic approach is to consider the linear member to be of higher order 
but still predominantly second order. Figure 22 shows the block diagram 
of a third-order system of this nature. The linear feedforward member 
could be, for example, the servo of figure 1 including the effects of. 
armature inductance. The same second-order feedback switching circuit 
was still employed. However, it could not be expected that the Pm^7n 
parameter values remain the same. Figure 23 compares the response of a 
third-order linear and nonlinear system. Here it is seen that the non- 
linear system still responded with much less error than the linear system. 

On the basis of the comparative performance of these two systems it 
seems important to further studies toward control of general higher order 
systems using techniques similar to those developed in this investigation. 


CONC HIDING REMARKS 


From control equation ( 8 ) a second-order nonlinear control system 
that tends to maintain small instantaneous error for relatively arbitrary 
inputs has been synthesized using digital-computer techniques. The only- 
type of nonlinear operations required in the realization are switching- 
type operations (zero-coincidence detection and parameter switching). 

The switching requirements are severe from the aspect of detector sensi- 
tivity and switching time delay but not impractical since all the switching 
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is done in feedback paths at low power levels and may thus be performed 
electronically. 

The system demands sensing of error and error rate of change. Since 
general noise in a system of this type has to be expected, both error 
and error rate are smoothed. In an extension of this work imperfections 
such as a definite noise level, a given threshold or time delay in the 
relays, and an overall dead time are being investigated more intensively. 

It is expected that they do not seriously impair the working of the system. 

It is a difficult task to compare the performance of a nonlinear 
system with that of a linear system because no general criterion for com- 
parison is available. Since for nonlinear systems the law of superposition 
does not hold, it is not adequate to choose the response to a certain 
input (e.g., the step input) as a criterion for performance comparisons. 

A number of different inputs therefore have been chosen for exhibiting 
the qualities of the nonlinear system. 

Experimental results indicate that this type of nonlinear control 
system performs better than a linear control system having a normalized 
frequency 15 times greater. Performance is evaluated in terms of the 
average value of the magnitude of the instantaneous error for band-limited 
inputs. Further, the nonlinear system performance is virtually independent 
of variations in the damping factor of the system. 


Stanford University, 

Stanford, Calif., November 23, 1955* 
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APPENDIX A 

COMPARISON OF EQUATION (8) AND NONLINEAR CONTROL EQUATION 
OF REFERENCES 1 AND 5 TO 5 

Comparison of Notations 

In original notation the nonlinear control equation was given els 

(refs. 1 and 3 to 5) 



(AI) 


where 


y 0 ( t ) 

yi(t) 


output 


input 


b = b + 1 + 


ZAJj -j 

— sgn (y^E) 

b 


Abp 

+ — — sgn 
b + 



X Ac-i . . 

C = c + 1 + sgn (y 0 E) 



E = (y Q - y±) 

( ) ' = d( )/dt 


Ab-j^AbgjAcipAcp constants, Ab-j_ > Abp, Ac-j_ > Acp 
a + ,b + ,c + constants 

In terms of the notation of equation (8) the above equation is written 



(A2) 
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where 

Pm = "lP s S n ^ y ' e ^ " 2P s S n (y' e ’)> m * 0. 1, 2, 5 

7 n = “l 7 sgn (ye) “ sgn (y e ')> n = 0, 1, 2, 3 


e = (x - y) 

pPjpP'i 7 ^ 7 

constants, > 2 P> 

l7> 2 7 

a + ,b + ,c + 

constants 


The subscript convention is 



P 3 - qP + 2 P 

73 = X 7 + 


P 2 = lP " 2 P 

OJ 

1 

K 
1 — 1 

II 

OJ 


h - 

7i = -7 2 


3 0 = -Pj 

7 0 = -73 


Evidently, the correlation between equations (Al) and (A2) is then 

y = y Q 

x = y ± 

Ah -1 Ah q Ac 1 Ac o 

1.0 - 2 0 = -r- ir - — » 2>~~r 

b b c c 

e = -E 

and a + , b + , and c + are the same. 

In either notation a set of seven (constant) parameters is needed 
to characterize the system. For example, in equation (A2) the set 

(a + , b + , c + , L P, 2 P, -J, 2 7) is sufficient. 
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Normalization of Nonlinear Control Equation 


When considering control systems it is possible to reduce the num- 
ber of parameters necessary to specify the nonlinear system. To do this 
a form of normalization familiar to linear theory is utilized. First, 
it should be realized that if error is to be defined els e = (x-y) 

then Cq = 1. Thus, assuming C + = 1 equation (A2) may be written 


Thus , knowing the natural frequency of the undamped linear system oo v , 

the number of parameters necessary to specify performance in the nonlinear 
case is five, that is. 



(A3) 


where 




Introducing normalized time r = o> v t, equation (A3) becomes 




(d, 1 s, 2 e, j?, 2 r) 


or alternatively 


(P> $2* ? 2 ) 
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APPENDIX /B 

SIMULATION OF SECOND-ORDER NONLINEAR CONTROL SYSTEM 


Equipment 

Experimental studies were carried out with the following equipment: 

( 1 ) General equipment : 

(a) Analog computer, Beckman Ease 

(b) Low-frequency function generator, hp 202 A (input device) 

(c) Pen recorder, Sanborn Twin-Viso (output device) 

(d) Vacuum-tube voltmeter, RCA WV-97A 

(e) Oscilloscope, Dumont 304-A 

(2) Special equipment: 

(a) Binary-logic switching circuit employed in conjunction 

with (a), described in detail in section "Binary-Logic 
Relay Switching Circuit." 


Computer Setup 

The computer diagram for the differential equation 

~ + 2D(1 + Pm) + (1 + 7 n )Y = *( T ) (Bl) 

dT 


is given in figure 24. The correspondence between this simulated system 
and the block diagram of figure 4 or 5 is straightforward (see ref. 6 ). 


Operational amplifiers 



through 



are used in simulating the 


linear member of the physical nonlinear system. Resistor R-j_ provides 

adjustment of the linear damping factor D. The input to this simulated 
linear member is x - (2 D(3 m y T + y n y) where x is obtained from the input 

device, while the values of feedback 2Dp m y' and y n y are obtained with 
resistors Rg through R^ (see table 5) connected to y T and y through 
a binary-logic relay switching circuit derived in the following section. 


Amplifiers 


through 


are sign changers. The four variables 


y, y ! , e, and e T whose signs are to be sensed are made available as 
shown in the lower right of figure 24. 
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Binary-Logic Relay Switching Circuit 

To complete the simulation of equation (8) it is necessary to imple- 
ment the binary logic of table 2. Figure 5 shows a relay switching cir- 
cuit (designed on an "and/or" basis) that realizes the necessary logic. 

In order to preclude the possibility of time delay in switching, 
fast" relays have been employed. The average properties of the 14 
double-pole single-throw relays comprising the switching circuit are 
given in table 4. 

The necessary synchronism in relays connected by dashed lines 
(fig. 5 ) is obtained by series connection of field coils. Into the four 
channels of field coils thus formed is read the four-digit binary logic 
of table 2 where now 

0 = Function > 0 = No coil current 

1 = Function < 0 = Coil current 

This process of reading in the binary logic may be done by monitoring 
the variable (y, y', e, e') with four separate zero-coincidence detectors 
or amplitude selectors (denoted CD in fig. 5 and described in detail 
in the following section). The output of these coincidence detectors 
then drives the respective channels of relay coils. 


Sensitive Transistorized Zero-Coincidence Detector 

One of the four identical amplitude selectors designed for use with 
the switching circuit of figure 5 is shown in figure 25. This circuit 
basically consists of a grounded emitter PNP transistor stage Tq driving a 

pentode pulse amplifier T 2 . (it should be possible to utilize a tran- 
sistor in place of the pentode if desired.) Four relay coils connected 
in series form the plate load of T 2 . Positive feedback R^ has been 
incorporated for regenerative switching. 

The operation of this circuit is as follows: When the input volt- 

age V becomes more negative than -50 millivolts, current starts to flow 
in the base-emitter region of Tq. This initiates collector-emitter cur- 
rent which effectively grounds the collector and thus the grid of T 2 . 

Plate current flows in T 2 so that the relays which form this plate load 
are switched to the up position (fig. 5) . Positive feedback Rj has been 
incorporated to make the switching regenerative. Stages Tq and T 2 will 
continue to conduct until V goes positive by 50 millivolts, at which 
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time current is cut off in T-^ and thus T2; the relays return to the nor- 
mally closed positions (down in fig. 5)* 

Refinements that have "been included in this circuit are: 

(1) Adjustable positive feedback R^. This feedback gives control 

over the zero-sensing threshold of the detector. The greater the posi- 
tive feedback the greater the threshold. This adjustment is desirable 
for studying the effects of threshold in sensing discussed in the section 
’’Effects of Switching Imperfections.” The measured peak-to-peak thresh- 
old values obtainable with this arrangement were a maximum of kOO milli- 
volts and a minimum of 60 millivolts, (in the computer setup it was 
necessary to amplify error e that formed the input to one of these 
detectors five times in order to bring the switching threshold down to 
15 millivolts, peak to peak.) 

(2) Adjustable bias for emitter of T-^. This is necessary to com- 
pensate for the slight positive bias (approximately equal to 100 milli- 
volts) given to the base of T^ by the positive feedback R^. 

(3) A clamping or clipping diode T^ to protect the transistor T^ 

from excessive base-emitter inverse voltages. Thus it can be seen that 
the input Impedance of this coincidence detector is 15 kilohms since the 
base of T^ is effectively always grounded. 
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APPENDIX C 

BEHAVIOR OF OUTPUT IF INPUT AND INPUT DERIVATIVE 
ARE SIMULTANEOUSLY SMALL 


Assume that the input is a curve 'which may he approximated by 

x = OCT ^ 

for values t <5. Since x* = 2ar, both input and input derivative 
vanish with t — > 0. The input phase curve is 

X = Mx') 2 

4a 

The output depends on the initial conditions: 

(1) For y(0) = 0 and y'(0) = 0 

t 4 D(l + ^ 

y = a — + — T'’ + 

J \12 30 

(2) For y(0) = 0 and y'(0) = e-^ 

y = - D(1 + pjr 2 + . 

(3) For y(o) = e 2 and y'(0) = 0 

y = e, 




1 + I 


1 T 2 + . . . 


V 2 / 



In all cases reproduction of the input is not perfect very close to 
t = 0; however, in cases (2) and (3) switch points may occur for rather 
small values of t. In case (l) no switch point close to t = 0 can 
be expected. This case will rarely occur; in most cases neither y(0) 
nor y'(0) will be zero. Then the output is a superposition of cases (2) 
and ( 3 ) • In this event the error e = x - y is given by 


e = -e 2 - + 


a + 6 ]_D(l + ^ - e 2(— + * * • 
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It is evident that for small values of e^ and eg the error grows 

with a. This can easily he seen in the results of figure 11. The input 
was x = A( 1 - sin At); that means , near x = 0, the input may he approxi- 
mated hy a parabola with a = l/2(ft 2 )A or the errors near x =* 0 grow 
with ft 2 . 
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APPENDIX D 
APPROXIMATING CURVES 


The approximating curves which can be used to trace the output for 
a given input form a network in the phase plane. It has proved practical 
to present the network for e and e 1 of equal sign in one sheet 
( e -g •) see fig. 20) and the network for e and e' with opposite sign 
on another sheet. Superposition of both sheets allows establishing of 
the four approximating curves through each point of the phase plane. 

The approximating phase curves are determined by equation (2l). 

Since the values and A 2 change from quadrant to quadrant the curves 

are composed of portions of different analytical curves which are patched 
at y' = 0 and y = 0. 

The roots A-^ and A 2 depend on D, p m , and 7 n j they may be 

complex or real. If the roots are complex, the approximating curves are 
curves of the spiral type which wind around the origin of the phase plane. 
In the case of real roots the curves have a quite different character. 

This can easily be understood by transforming the equation of the approxi- 
mating curves (see refs. 1 and 3 to 3 and ch. V of ref. 7). For real 
roots the new coordinates p and v are introduced: 

yAl - Ip = p 

y\ 2 _ I - ! = v 

Then equation (2l) yields 



If A^ and Ag are real and of opposite sign, the curves in the 

pv plane have hyperbolic character with (p,v) = (0,0) as saddle point 
(see fig. 26(a)). If Aj_ and A2 are of equal sign, the curves have 

a nodal point (see figs. 26(b) and 26(c)). In figures 27 and 28 phase 
curves with saddle point and stable and unstable node are shown in the 
original yf^ plane. 

The approximating curves are composed of portions of these different 
types (see fig. 29). 
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There is no need of avoiding |3 ni ,7 n combinations which lead to 

node-type approximating curves because only portions of these curves 
are used. In the earlier publications it appeared as if (in either one 
of the networks of approximating curves) one set of curves had to be 
formed by curves of spiral character (complex roots A). However, this 
has proved to be an unnecessary restriction. There might be some trouble 
with node- type curves if large delays in switching should occur (e.g., 
delays in y' = 0 switching in fig. 29). 
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TABLE 1.- CODING SCHEME FOR SWITCHING LOGIC FORM EQUATION (l4) 
[b, function > 0; 1, function < 6] 


Binary coded decimal 

Decimal 

y 

y' 

e 

e 1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

0 

0 

1 

0 

2 

0 

0 

1 

1 

3 

0 

1 

0 

0 

k 

0 

1 

0 

1 

5 

0 

1 

1 

0 

6 

0 

1 

1 

1 

7 

1 

0 

0 

0 

8 

1 

0 

0 

1 

9 

1 

0 

1 

0 

10 

1 

0 

1 

1 

11 

1 

1 

0 

0 

12 

1 

1 

0 

1 

13 

1 

1 

1 

0 

14 

1 

1 

1 

1 

15 
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TABLE 2.- MATRIX OF ALLOWED ^ m ,7 n COMBINATIONS AS DETERMINED 
BY ENCODED SWITCHING LOGIC OBTAINED FROM EQUATION (lL) a 



3o 

Pi 

02 

P 3 

?0 

0000 

1111 



0100 

1011 

7l 


0001 

1110 

0101 

1010 


72 


0110 

1001 

0010 

1101 


73 

0111 

1000 



0011 

1100 


a ExampIes: Pj7q is chosen when 0100 

occurs, i.e., 

(y > 0, y 1 c 0, e > 0, e' > 0) 
or when 1011 occurs, i.e., 

(y < 0, y’ > 0 , e < 0 , e' < 0) 

P 2 7i is chosen when 0101 occurs, i.e., 

(y > 0, y' < 0, e > 0, e' < 0) 
or when 1010 occurs, i.e., 

(y < 0, y' > 0, e < 0, e' > 0) 


is not possible. 
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TABLE 3-- PARAMETERS FOR FIGURE I 6 


Figs. 16(a) and 16(b) 

Figs. l 6 (c) and 16 (d) 

Complete system 

System 1: 

P3 = -P 0 =2 7 ^ = - 7q = 2 

p 2 = = 0.5 7 2 = -7 L = 0-5 

System 5*- 

P3 = “Pq = PO 73 = ~7q = 10 

P 2 = -Pp =1 72 = "7i = 1 

(l) No derivative feedback 

System 2: 

P3 = -Pq =0 = - 7 q = 2 

= "Pi = 0 72 = “ 7 l = 0-5 

System 6 : 

P3 = -Pq =° 73 = “7q = 10 

p 2 = -p L = 0 7 2 = - 7 X = 1 

(2) No e' sensing 

System 3: 

P 3 = -Po = 2 73 = -7o = 2 

P2 = “Pi = 2 7 2 = “7p = 2 

System 7 : 

P3 = “Po = 10 73 = - 7 o = 10 

?2 — “Pp = 10 72 = “7 1 = 10 

(3) No e sensing in derivative feedback loop and 
no e' sensing in proportional feedback loop 

System k: 

P3 = -Po = 2 73 = -7q = 2 

P 2 = -p ] _ = -2 7 2 = “7 ]_ = 2 

System 8 : 

Pj = "Pq = ^ = ~ 7 q - 10 

P2 = “Pp = ~10 72 = ~7p ~ 10 
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TABLE 4.- AVERAGE PROPERTIES OF RELAYS FOR SWITCHING 
CIRCUIT OF FIGURE 5 


Description: General Electric CR-2791 

double-pole single-throw relay mounted 
on 5~P r °ng Amphenol base with permanent 
aluminum dust cover 


Connections: 



Electrical Properties 

Values 

Coil resistance 

2-5 kilohms 

Coil current for 


positive action 

5 to 10 ma 

Pull-in time 

2 millisec 

Drop-out time 

3 millisec 
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TABLE 5.- CORRESPONDENCE BETWEEN ADJUSTABLE COMPUTER 
ELEMENTS AND PARAMETERS OF EQUATION (8) 


Element 

Corresponds to 

Units , 
kilohms 

R i 

( 2 D) -1 

100 

r 2 

(DDPq)- 1 

100 

R 3 

( 2 DPJ- 1 

LOO 

r 4 

(2DS 2 )-1 

100 

R 5 


100 

r 6 

1 1 

1 

O 

100 

*7 

r — 1 

1 

1— 1 

100 

r 8 

1 — 1 

1 

OJ 

100 

*9 

(73)- 1 

100 
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Figure 1.- Block diagram of simple positional servo. J, inertia of motor 
rotor, gears, and reflected load; B, viscous damping of motor and 
reflected load referred to motor shaft; k^, constant of proportionality 

between output velocity and back electromotive force (it includes arma- 
ture resistance and gear ratio from motor shaft to load). 



7 


n 




>7 


Gain 


7- 


7c 


7l 


- 7 °J- 


Figure 2.- Illustration of stepwise nature of parameters 3m and / n 
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Figure 5 



\ ^ 


- Block diagram of equation (8) assuming that 0m and 
constant (denoted by encircling dotted line). 


7n 


are 
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logic | 

1 1 


Figure 4.- Block diagram of equation (8) taking into account stepwise 

switching nature of and j n . 
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y 

Output 


Figure 5 .- Block diagram showing complete physical interpretation of 
equation (8). CD, zero-coincidence detectors. 
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(a) Nonlinear system; ft = od/cjd v = 0.1 to 0.4 in 0.1 steps. 



(b) Linear system; ft = cd/o^ = 0.1 to 0.4 in 0.1 steps. 

Figure 6 .- Linear and nonlinear system responses for 20-volt peak-to- 
peak sinusoidal inputs with frequency ft varied. D = 0.6; 2.5 small 
divisions on time scale = 1 normalized time unit; tick marks at 
bottom of figures indicate where frequency was varied. 
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(c) Nonlinear system; fl = (u/auy = 0.5 to 0.8 in 0.1 steps. 



(d) Linear system; 0 , = oo/ov = 0.5 to 0.8 in 0.1 steps. 
Figure 6.- Continued. 
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(f ) Linear system; ft = o>/a>v = 1.0 to 2.0 in 0.2 steps 
Figure 6.- Concluded. 


(a) Nonlinear system; £1 = 2tt/T = 0.1 to 0.4 in 0.1 steps. 



(b) Linear system; 0 , = 2jc/T = 0.1 to 0.4 in 0.1 steps. 

Figure 7 .- Linear and nonlinear system responses to 20-volt peak-to-peak 
triangular -wave inputs with period T varied. D = 0.6; 2.5 small 
divisions on time scale = 1 normalized time unit; tick marks at bottom 
of figures indicate when periods were varied. 


NACA TN 3826 


61 



(c) Nonlinear system; 0 . = 2rt/T = 0.5 to 0.8 in 0.1 steps. 



(d) Linear system; ft = 2n/T = 0.5 to 0.8 in 0.1 steps. 


Figure 7 .- Continued 
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(f) Nonlinear system) Q = 2jt/T = 1.0 to 2.0 in 0.2 steps. 


Figure 7*- Concluded. 
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(a) Nonlinear system; a = 4 and 2. Tick mark at bottom of figure indicates 

where a was varied. 



(b) Linear system; a = Ij- and 2. Tick mark at bottom of figure indicates 

where a was varied. 

Figure 8.- Linear and nonlinear system responses to 20-volt peak-to-peak 
partially integrated square-wave input. D = 0.6; 2.5 small divisions 
on time scale = 1 normalized time unit. 
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(c) Nonlinear system; a = 1 and 0.5. Tick mark at bottom of figure 
indicates when a was varied. 



(d) Linear system; a = 1 and O. 5 . Tick mark at bottom of figure indicates 

when a was varied. 


Figure 8 .- Continued 
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(e) Nonlinear and linear systems; a = 0. Tick mark at bottom of figure 
indicates demarcation between nonlinear and linear systems. 


Figure 8.- Concluded. 
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(b) Linear system; = 0.1 to 0.4 in 0.1 steps. 

Figure 9-- Linear and nonlinear system responses to 4-volt peak-to-peak 
sinusoidal inputs with frequency fi varied. D = 0.6; 2.5 small 
divisions on time scale = 1 normalized time unit; tick marks at 
bottom of figures indicate where frequency was varied. 
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(d) Linear system; SI = o)/co v = 0.5 to 0.8 in 0.1 steps. 


Figure 9.- Continued 
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(e) Nonlinear system; fi = cd/o^ = 1.0 to 2.0 in 0.2 steps. 



(f) Linear system; 0, = cojux^ = 1.0 to 2.0 in 0.2 steps, 
Figure 9-- Concluded. 
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(a) Nonlinear system; ft = a = 0.1 to 0.4 in 0.1 steps. 



(b) Linear system; ft _ co^GUy = 0.1 to 0.4 in 0.1 steps. 

Figure 10.- Linear and nonlinear system responses to 20-volt peak-to- 
peak sinusoidal inputs with frequency ft varied that are clipped 
symmetrically at ±6 volts. D = 0.6; 2.5 small divisions on time 
scale = 1 normalized time unit; tick marks at bottom of figures indi- 
cate where frequency was varied. 
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(d) Linear system; ft = atjox^ = 0.5 to 0.8 in 0.1 steps. 


Figure 10.- Continued 


MCA TN 3826 


71 



(f) Linear system; 



1.0 to 2.0 in 0.2 steps. 


Figure 10 


Concluded 



72 


NACA TN 3826 




Figure 11.- Nonlinear system responses of 20-volt peak-to-peak sinusoidal 
input with frequency Cl varied that has been displaced by -10 volts 
(direct -current component) for obtaining osculation. D = 0.6; 2.5 small 
divisions on time scale = 1 normalized time unit; tick marks at bottom 
of figures indicate where frequency was varied. 
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(b) Linear system. 

Figure 12.- Linear and nonlinear system responses for triangular-wave 
input whose period and amplitude are randomly modulated. D = 0.6; 
1 small division on time scale = 1 normalized time unit. 
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(a) Nonlinear system response to triangular -wave input with period varied 
as in figure 7(a) • D = 0.6; tick marks at bottom of figure indicate 
when period was varied. Threshold in sensing sgn e was 44 millivolts. 



(b) Constant input x(t) was 9 volts. 


Figure 13*- Effects of switching delays due to threshold in sensing sgn e. 
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(b) Nonlinear system; ±12-volt y" limit. 

Figure 14.- Nonlinear system responses to 20-volt peak-to-peak sinusoidal 
input with varying acceleration limits. D = 0.6; SI = o)/go v = 0.5; 

2.5 small divisions on time scale = 1 normalized time unit. 
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(a) 


Nonlinear system; ±7 


1 

2 


-volt 


limit . 


Figure 14.- Continued 
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(e) Nonlinear system; ±6-volt y" limit. 



(f) Nonlinear system; 4|-volt y" limit. 


Figure l4 


Concluded 
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(b) Linear system. 

Figure 15.- Linear and nonlinear system responses to 20- volt peak-to- 
peak sinusoidal input with ±4^ -volt velocity limit. D = 0.6; 

ft = 2rt/T = 0.5; 2.5 small divisions on time scale = 1 normalized 
time unit. 
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(a) Complete system 1 and system with no derivative feedback 2. 

System 1: ^ = -p Q = 2; p 2 = - Pl = 0.5; 7? = - 7q = 2; y Q = - 7± = 0.5. 

System 2: ^ = -p Q = 0; P 2 = = 0; 7? = -y Q = 2; 7£ = - 7± = 0.5. 



(b) System with no e’ sensing 3 and system with no e sensing in deriva- 
tive feedback loop and no e f sensing in proportional feedback loop k. 
System 3 : „ } = . 2; % = - 3l - 2; . -y Q . 2, y s . , 2 . 


System 


h m 


■P 0 “ 2 > P 2 ” _(3 1 " ~ 2 ' ^ = - 7q = 2; y = - 7 


= 2 . 


Figure 16.- Responses of nonlinear system with four different parameter 
sets to a 20-volt peak-to-peak triangular-wave input with fixed period 
10*. D 7 °‘°i 2 *5 small divisions on time scale = 1 normalized 
time unit; tick marks at bottom of figures indicate demarcation between 
systems • 
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(c) Complete system 5 and- system with no derivative feedback 6. 

System 5: (3^ = -0 Q = 10; 0 2 = = 1; y^ = -y Q = 10; j 2 = -y^ = 1. 

System 6: (3^ = -0 Q = 0; 0 2 = = 0; y^ = -7 Q = 10; y 2 = -7 X = 1. 



(d) System with no e' sensing 7 and system with no 
tive feedback loop and no 0 ' 

System 7 : 

System 8: 


e sensing in deriva- 


- -U 0 = 10! 


f>3 - -l»o ‘ 10; ‘ 


e 1 
>2 
3 2 


sensing in proportional feedback loop 8. 
x = 10; 7^ = -7 q = 10; 7 g = -y ± 

3 _ _/ 0 


0 X = -10; 7* = -7 0 = 10 ; 7 2 = "7 1 = 10 • 


Figure l6.- Concluded 
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(a) fi = 2n/T = 0.1 "bo 0.4 in 0.1 steps. 



(b) Q, 2ji/t — 0.5 to 0.8 in 0.1 steps. 

Figure 17.- Nonlinear system responses to 20-volt peak-to-peak triangular- 
wave input with period T varied. D = 0.1; = _p Q = 12; p 2 = _p = 


7* = -7, 


0 


2 » 7 2 7j_ - 0.3; 2.5 small divisions on time scale = 1 nor 

was Varied 1116 markS at bottom of figures indicate where perio 
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Figure l8.- Sketch showing portions of superposed input and output of 
nonlinear system magnified approximately 100 times. 
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Figure 19.- Phase-plane angular sectors defined by equation (22) for the 
given parameter values. 1 P = 1.5; = 0.5; = 0.^5; - y = 0.05; 

2D = 0.25; % = 2. 



Figure 20.- Superposed input and output phase planes showing available 
angular sectors for sinusoidal inputs, ^p = 1 . 5 ; ^7 = 0.45; ^P = 0 . 5 ; 


2 y = °* 0 ^ 20 = °’ 2 ^ = 2 > Mmax = 2 ' Mma* = 


I min 


= i; 


71 . =0.4. 

' Imin 
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Figure 21.- Portion of output phase diagram for square -wave input, 

P 3 = _(3 0 = 2 ’ = " P 1 = °'^> = _7 o = 2 > ?2 = " 7 1 = 0, 5> D = °*6j 

o\ = 1 . 
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Figure 22.- Block diagram of third-order nonlinear control system. 
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(b) Linear system. 

Figure 23-- Linear and nonlinear third-order system responses to random 

input. T = 0.1; D = 0.6; 2.5 small divisions on time scale = 1 nor- 
mail zed time unit. 
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Figure 2k. - Computer diagram for simulation of equation (8). 
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Figure 25-- Transistorized zero-coincidence detector. 
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Figure 26.- Curves in |av plane. 
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Figure 28.- Phase curves with nodal point. 
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Figure 29.- Composition of an approximating curve. 
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